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We introduce a set of coherent states which are associated with quantum 
systems governed by a trilinear boson Hamiltonian. These states are produced 
by the action of a nonunitary displacement operator on a reference state and 
\q . can be equivalently defined by some eigenvalue equations. The system pre- 

pared initially in the reference state will evolve into the coherent state during 
the first instants of the interaction process. Some properties of the coherent 
states are discussed. In particular, the resolution of the identity is derived and 
the related analytic representation in the complex plane is developed. It is 
, shown that this analytic representation coincides with a double representation 

based on the Glauber coherent states of the pump mode and on the SU(1,1) 
Perelomov coherent states of the signal-idler system. Entanglement between 
the field modes and photon statistics of the coherent states are studied. Con- 
nections between the coherent states and the long-time evolution induced by 
. the trilinear Hamiltonian are considered. 
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I. INTRODUCTION 
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Various nonlinear processes in quantum optics, such as parametric amplification, fre- 
quency conversion, Raman and Brillouin scattering, and the interaction of two-level atoms 
with a single-mode radiation field, are described by a trilinear boson Hamiltonian of the 
form jg-gJJ 



H = uj a b)a + u$b + u c c^c + koP$ + tCa^bc. (1-1) 



Here a, b, c are the boson annihilation operators of the three field modes with angular 
frequencies u a , cj&, u c , respectively, which satisfy the energy conservation law 

u a = u b + u c . (1.2) 



The concept of coherent states has been introduced by Glauber |3(J for the single-mode 



quantized radiation field whose dynamical symmetry group is the Weyl-Heisenberg group 
f31"| . Perelomov |32|] and Gilmore f33| , |34|j have generalized the concept of coherent states to 



quantum systems whose dynamical symmetry group is an arbitrary Lie group. 

In the present paper we introduce coherent states for quantum systems whose dynamics 
is governed by the trilinear boson Hamiltonian fll.ip . This Hamiltonian is a linear combi- 
nation of operators that do not form a finite-dimensional Lie algebra, and hence a direct 
application of the Perelomov formalism becomes impossible in this case. Therefore, we de- 
fine the coherent states as eigenstates of some non-Hermitian operators or, equivalently, as 
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states produced by the action of a nonunitary displacement operator on a reference state. 
We show that during the first instants of the time evolution (i.e., in the short-time approxi- 
mation) this reference state is evolved into the coherent state. Moreover, the coherent states 
have an important property of resolution of the identity. Thus we are able to develop an 
analytic representation in the complex plane. Some quantum statistical properties of the 
coherent states are also discussed. In particular, we study photon statistics of the pump 
mode and the entanglement between the pump and the signal-idler subsystems. We also 
consider connections between the coherent states and the long-time dynamics governed by 
the trilinear boson Hamiltonian. 



II. PHYSICAL MODELS 

The dynamics induced by the Hamiltonian ( |1.1| ) can be better understood by revealing 
its symmetries. It can be easily shown that the operators 

S ab = a} a + Pb, (2.1a) 
Sac = a^a + c^c, (2-lb) 
V hc = Pb-$c (2.1c) 

(known in literature as Manley-Rowe invariants) are integrals of motion. They manifest 
underlying SU(2) and SU(1,1) symmetries. Only two of these invariants are linearly inde- 
pendent: 

S a b ~ S ac = T> bc - 

Therefore, a state of the system can be characterized by two quantum numbers. The Hamil- 
tonian (II . II) can be written in the form 



H = H + H int (2.2) 

where 

H = uj a a)a + tUbWb + cu c c^c, (2.3a) 

H int = K aPc ] + n*a r bc, (2.3b) 

with 

[H ,H iQt ] = 0. (2.4) 

With the choice u b = u c = u a /2, the Hamiltonian ( |I~T| ) can be written (up to an 



unimportant constant) as 

H = uj a (a)a + K ) + KaK + + /t*a f K__, (2.5) 

where the operators 

K + = Pc\ = be, K = \{Pb + c f c + 1) (2.6) 
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span the SU(1,1) Lie algebra: 



[k , k ± ] = ±k ± , [k_, k + ] = 2k . (2.7) 

The form ( |2.5| ) is convenient for description of processes with underlying SU(1,1) symmetry, 
e.g., for parametric amplification where the modes d, b, and c are identified as the pump, 
signal, and idler modes, respectively. The signal-idler system possesses an SU(1,1) symmetry, 
and the Casimir operator 

k 2 = ki-\(k+k- + k-k+) = k{k-\) (2.8) 

is an integral of motion: 

K 2 = \T>l-\. (2.9) 

Thus the Bargmann index k that labels irreducible representations of SU(1,1) is related to 
the photon-number difference between the signal and idler modes: 

k = \{\V bc \ + 1). (2.10) 

Thus k takes discrete values k = |,2, ... that corresponds to the so-called discrete 
series of representations |3l| . The representation Hilbert space is spanned by the complete 
orthonormal basis \k,n) (n = 0, 1, 2, . . .), 

K \k,n) = (k + n)\k, n). 

For the two-mode realization ( |2.6j ) the states \k, n) can be expressed in terms of Fock states 
of the modes: 



\k,n) = \n + 2k-l) b \n) c . (2.11) 

This form implies that there is more photons in the mode b than in the mode c. However, 
since the Hamiltonian is symmetric in the operators b and c, we may by definition refer to 
the mode with larger number of photons as the mode b and to the mode with smaller number 
of photons as the mode c. Certainly, this definition is incomplete if an initial state is given, 
for example, by a superposition |n)&|0) c + |0)b|m) c . In such a case we must use two SU(1,1) 
representations: \k,n) + = \n + 2k — l)j,|n) c and \k,n)- = \n) b \n + 2k — l) c . However, for 
the sake of simplicity, we will not consider such situations. 

There is a number of nonlinear optical processes described by the Hamiltonian ( |2.5| ). For 
degenerate parametric amplification or second harmonic generation, one has 

k + = \P\ k_ = \b\ k Q = \$i> + \). (2.i2) 

In this case K 2 = —3/16, so k = 1/4 or 3/4. The dynamics induced by the corresponding 
Hamiltonian has been studied for both second harmonic generation p6| , |37| and degener- 
ate parametric amplification |38|-4l2]]. For the intensity-dependent two-mode coupling, the 
SU(1,1) generators are given by 
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k + = \Jm h\ k_ = hytib, k = tib + §, (2.13) 

which is known as the Holstein-Primakoff realization [4H,JI|. In this case k = 1/2. 



The Hilbert space Ti of a quantum system governed by the Hamiltonian (|2.5| ) can be 
decomposed into a direct sum of finite-dimensional subspaces TCl- 

oo 

H = @H L . (2.14) 

L=0 

Each subspace TCl is spanned by the complete orthonormal set of (L + 1) states \n)\k, L — n) 
(n — 0, . . . , L), where \n) are the Fock states of the mode a, 

a*a\n) = n\n), 

and \k,L — n) are the orthonormal states of SU(1,1). We have seen that the Bargmann 
index k is related to the integral of motion Z>& c via Eq. (|2.10|) . Correspondingly, we find the 
relation between the integrals of motion and the quantum number L: 

L = \{S ab + S ac -V bc ). (2.15) 

It is clear that the operator 

H = Lu a (tia + K ) (2.16) 

has the property 

H \n)\k,L-n) — u a (k + L)\n)\k, L — n). (2.17) 
Since Hq commutes with 

H int = KdK + + K*tfk_, (2.18) 
the subspace TCl is invariant under the evolution induced by the Hamiltonian fl2~5|) . 

III. COHERENT STATES 

The evolution operator U(t) = exp(-iHt) can be written as 

U(t) = exp(zaK + - z*a ] K^) exp(-iH t), (3.1) 

where z = —int. If the initial state is \ipQ\ k, L), the action of exp(— iH^t) will only multiply 
it by a phase factor: 

exp(-iH t)\ip ; k,L) = exp[-iu a (k + L)t}\ip ; k,L). 
Therefore we can consider only the action of the first exponent in ( |3.1|) : 

\ip z ;k,L) = exp(zaK + - z*a [ kJ)\il)Q\ k, L). (3.2) 



4 



If the mode a is highly excited and its depletion can be neglected through the evolution of 
the system, one can use the so-called parametric approximation I n this approximation 

the mode a is treated classically, i.e., the operators a and a) are replaced by c numbers a 
and a*, respectively. In this case Eq. ( |3.2|) defines the SU(1,1) Perelomov coherent states. In 
the context of the boson realizations ( |2.6| ) and ( 2.12j ), these states are identified as the well- 
known two-mode and one-mode squeezed states, respectively. (Limitations to squeezing 
in degenerate parametric amplifiers due to quantum fluctuations in the pump have been 
discussed in Refs. pS-fjL?]].) The SU(1,1) coherent states have been also considered [flfi-l? 
in the context of the Holstein-Primakoff realization ( |2.13| ). 

In what follows we will consider the general (pure quantum) case, when the parametric 
approximation is not valid. However, it is customary to use the short-time approximation 
p!0| , p!2|JT7| , |2^p7| , |29[] in which one assumes \z\ <C 1. In order to clarify the algebraic mean- 
ing of this approximation, we use the Baker- Campbell-Hausdorff theorem expressed as the 
Zassenhaus formula [ ff8| : 

Here Cj is a polynomial of jth degree in X and Y with rational coefficients. For example, 

C 2 = -±[X,f], (3.4) 
C, = -\[[X,Y],Y]-\[[X,nX]. (3.5) 



For the evolution operator of Eq. (3^), we identify 

X = zaK + , Y = -z*a)K_. (3.6) 
Therefore, Cj contains \z\K For example, 

I |2 I 12 

C 2 = ^[ak+itfk..] = ¥±-[k-K+ - 2atik }. 

In the short-time approximation one truncates in the second order of \z\, i.e., neglects all 
the terms that contain \z\ J with j > 2. In this approximation we obtain 

\ip z ; k, L) w [l + (zaK + - z*a^K_) + \(zaK + - z*a^k^) 2 ] \tfj ; k, L). (3.7) 
Let us consider the initial state of the form 



\fak,L) = \L)\k,0). (3.8) 



Then Eq. (|3.7|) gives 



\ip z ;k,L) w (1 - kL\z\ 2 )\L)\k, 0) + zV2kL\L - l)\k,l) 

+z 2 J(2k 2 + k){L 2 -L)\L- 2)\k, 2). (3.9) 



The initial state (|3.^ ) is an extreme state in the subspace Hl and hence it has some special 
properties: 



5 



e d2 \L)\k^) w (l + \\z\ 2 [aK + ,a)k^ \L)\k,0) = (1 - kL\z\ 2 )\L)\k, 0), (3.10) 



Therefore, we find 



exp(-2*a t i^_)|L)|A;,0) = |L)|fc,0). 



\if> z ]k,L) ss (1- A;L|2| 2 )exp(^ai' + )|L)|A;, 0), 



(3.11) 



(3.12) 



where the exponential should be expanded up to the second order, which immediately leads 
to Eq. (p|). 

Now we introduce the coherent states for the quantum systems governed by the Hamil- 
tonian (|2.5|) . We define 



\z;k,L) = e*e f e d2 \L)\k, 0) = 7V" 1/2 exp(zaK + )\L)\k, 0), 
where M is the normalization factor determined by 

exp(^\z\ 2 [ak + ,d jf k^}) |£)|Jfe,0) = N~ 1/2 \L)\k, 0). 



(3.13) 



Comparing Eqs. ( |3.13|) and ( 3.12 ), we see that in the short-time approximation the coherent 
states \z;k,L) coincide with the states \ip z ;k,L) obtained by the action of the evolution 
operator on the initial state \L)\k, 0). Expanding the exponential in (|3.13|) , we obtain 



z; k,L) = —= d n (k 1 L)z n \L - n)\k,n), 



(3.14) 



d n (k,L) 



LlT{2k + n) 



nl/2 



(3.15) 



_n!(L-n)!r(2A;) 

Sometimes it is convenient to use y = 1/z. Then we write the coherent states in the form 

-L L 



y; k,L) = ^ d n (k, L)y n \n)\k, L — n), 



(3.16) 



d n (k,L) = d L - n (k,L) 



L\T{2k + L-n) 
n\(L-n)\V(2k) 



n 1/2 



The normalization factor is given by 



Af(\y\ 2 ;k,L) 



T(2k + L) 
\y\ 2L T(2k) 



$(-m-L-2k-\y\ 2 ), 



(3.17) 



where $(a; c; x) is the confluent hypergeometric function (the Kummer function) 
scalar product of two coherent states is 



. The 
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{yi;k,L\y 2 ;k,L} - 



myi\ 2 ;k,L)Af(\y 2 \*;k,L)} 1/2 ' 

Equation ( |3 . 1 3| ) shows that the coherent states \z\ k, L) are produced by the action of the 
(nonunitary) displacement operator exp(zaK + ) on the reference state \L)\k,0). They can 
be equivalently defined by some eigenvalue equations. We start with eigenvalue equations 
satisfied by the reference state |L)|fc, 0): 

tfa\L)\k,0) = L\L)\k,0), (3.18a) 
k \L)\k,0) = k\L)\k,0), (3.18b) 
K_\L)\k,Q) = 0. (3.18c) 

Acting on both sides of these equations with exp(zak + ), we obtain the eigenvalue equations 
satisfied by the coherent states: 

(a)a + zaK + )\z; k, L) = L\z\ k, L), (3.19a) 
(k -zak+)\z;k,L) =k\z;k,L), (3.19b) 
(AT_ - 2zaK + z 2 a 2 K + )\z; k, L) = 0. (3.19c) 



Note that Eqs. (gl9aD and ( pISED are connected via the relation a) a + Kq = k + L. 



It is straightforward to generalize the concept of the coherent states to the case when 
the reference state of the pump mode is a superposition of the Fock states \L) with different 
values of L and the reference state of the signal-idler system is a superposition of the SU(1,1) 
states \k,0) with different values of k: 

oo 

|Vo> = £ W£<felM>- 

L=0 k 

We assume that |^o) is normalized. Then the coherent state is given by 

1 = exp(zaK + )|^o) = 

k,L 



\z) = -^=exp(zak + )\ip ) = ^g k h L \z;k J L). (3.20) 



For example, one can consider a case when the reference state of the pump mode is the 
Glauber coherent state or the squeezed state. In what follows we will study the coherent 
states with specific values of k and L. Properties of the coherent superposition (|3.20|) depend 
on the choice of the coefficients g k and hi. Expectation values for the state \z) are given by 
weighted sums over the corresponding expectations for the states \z\ k, L). 

IV. RESOLUTION OF THE IDENTITY AND ANALYTIC REPRESENTATIONS 

An important property of coherent states is the resolution of the identity. For the 
coherent states \y; k, L) we have 



J dfx{y;k,L)\y;k,L)(y;k,L\ = 1 L , (4.1) 
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where 1^ is the projection operator on the subspace Hl- 



n=0 



(4.2) 



In order to find the integration measure d[x(y; k, L) that realizes the resolution of the identity 
( |OD , we substitute y = y/re 1 ^, i.e., 

r = \y\ 2 , <P = argy, 

and define 

d/i(y; k, L) = R{\y\ 2 ; k, L)d 2 y = \R(r; k, L)drd<p. (4.3) 
The integration is over the whole complex plane. Integrating over the angle 0, we obtain 



l T(2k + L-n) 
h n\{L-n)\ 



r n T(r;k,L)dr 



where we have defined 



T{r;k,L) = 



ttL\ 



\n)\k, L — n) (n\{k, L — n\ = 1l, 



R(r;k,L) 



(4.4) 



T(2k + L) $(-£;! - L - 2k; r)' 



(4.5) 



By using the integral [^9], p. 285], 



o 



T(a)T(c — b) 



(4.6) 



we obtain 



T(r; k, L) = ^/^t^ : &(L + 2; 2k + L + 1; -r) 



T(2fc + L + 1 

This gives us the final expression for the integration measure: 
1 L + l 



dn(y; k, L) 



2nL + 2k 



$(-L; 1 — L — 2k; r)$(L + 2; 2k + L + 1; -r)drc 



(4.7) 



The resolution of the identity is important because it allows the use of the coherent states 
as a basis in the state space. Let us consider a state |/) in Hl'- 



n=0 



\f) = Y,fn\n)\KL-n). 
We introduce the analytic function 

Y f (y; k, L) = M 1/2 y L (f\y; ^> = E d n (k, L)f\ 



(4.8) 



L 

E 

?1=0 



nil i 



(4.9) 
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that determines the representation in the coherent-state basis: 



\f) = J df,(y-,k,L) Y ^^ \y-,k,L), (4.10) 



and the scalar product of two states is 



, , 4 r , ^YJy;k,L)Yt(y;k,L) 

(g\f) = / dn(y;k,L) ' ,) ■ ( 4 - n ) 



M\y\ 2L 

We also find the reproducing kernel of this representation: 

Y f (y; k,L) = J d»(x; k, L ) ^^k^j Yf(x; *' L) ' ^ 

The operators aK + , a^K-, and Kq act in the Hilbert space of analytic functions Y(y; k,L) 
as linear differential operators: 

aK + = -y£; + (L + 2k-l)± (4.13a) 

^K^ = -y 2 ± + L y, (4.13b) 
dy 

K = -y^- + k + L } (4.13c) 
dy 

and the expression for a)d follows immediately from the relation d^d + K = k + L. 
The coherent state \y ; k, L) is represented by the analytic function 

y L N(y* y;k,L) 

YM) -^(\y \^L)- (414) 

By using Eq. (|3.19a|) , we find that this function satisfies the second-order differential equation 
d 2 Y dY 

y V " [vlv +{L+2k ~ l)] ^ + Ly *° Y ° = °- (4 - l5) 



Eigenstates of the Hamiltonian (|2.5| ) are determined by 

H \v,k,L) =uj a (k + L)\u,k,L), (4.16) 
H int \v,k,L) = v\i/,k,L). (4.17) 

The corresponding analytic function 

Y v (y; k, L) = Af 1/2 y L (", k, L\y- k, L) (4.18) 

satisfies the second-order differential equation 

d 2 Y dY 
V^i + [iy 2 -(L + 2k-l)]-^ + (D- 7 Ly)Y v = 0, (4.19) 
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where 



7 = K/K , V = U / K . 



(4.20) 



The above results for the analytic representation in y plane can be easily transformed 
for the analytic representation in z plane. We just use the expression (|3.14j) for the coherent 
states and obtain 



Z f (z;k,L)=Ar^ 2 (f\z;k,L) = £ d n (k , L) f* L _ n z n . 

n=0 

It is easy to find the following relation: 

Z f (z;k,L) = z L Y f (y=l/z;k J L). 
The corresponding differential operators are 



aK 4 



d^K_ = - 

dz 

T> d , 

K = z— + k. 

dz 



-z 3 -^ + (L — 2k — l)z 2 4- + 2kLz, 
dz 1 dz 

d 



(4.21) 

(4.22) 

(4.23a) 
(4.23b) 
(4.23c) 



It is interesting to study how the coherent states \y; k, L) of the whole three-mode system 
are related to the Glauber coherent states of the pump and to the SU(1,1) Perelomov coher- 
ent states of the signal-idler system. We consider this relation by using the corresponding 
analytic representations. The Glauber coherent states are given by Q 



\a) 



co n 
| a |2/2^ a 



E 

n=0 



(4.24) 



and the corresponding analytic representation in a plane (called the Bargmann represen- 
tation) is well known |50|-p2|. The SU(1,1) Perelomov coherent states can be expressed as 



|fc,C) = (l-|C| 2 ) fc E 



n=0 L 



Y(n + 2k) 
n^(2k) 



1/2 



C\k,n), 



(4.25) 



where \(\ < 1. The corresponding analytic representation is referred to as the representation 
in the unit disk. We combine the Bargmann representation and the representation in the 
unit disk in order to obtain a double representation for the state |/) G TCl- 



D f (a,C) = e^\l ~ \t\r k (fm\k,0) = E [ X- + ^!r(2i) ] 1/2/ " aWC£ "" 
The differential operators of the double representation are given by 



(4.26) 
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d 

tfK- = a—, (4.27b) 
d( 

K = C-^ + k. (4.27c) 

It is easy to find the following relation: 

Df(a, C) = ^Y f (y = a/() = ^Z f (z = (/a). (4.28) 

This relation shows that the analytic representation in the coherent-state basis \y; k, L) is 

equivalent to the double representation based on the Glauber coherent states of the pump 
and on the SU(1,1) Perelomov coherent states of the signal-idler system. 



V. MODE ENTANGLEMENT AND QUANTUM STATISTICS OF THE 

COHERENT STATES 

The quantum nature of the dynamics induced by the trilinear boson Hamiltonian leads 
to strong entanglement between the pump and the signal-idler subsystems [ P7[| . We study 
here the entanglement between the field modes when the whole three-mode system is in the 
coherent state \y; k, L). This is a pure state and therefore the total entropy is zero. However, 
quantum correlations between the modes lead to the increase of the marginal entropy of each 
mode. The von Neumann quantum entropy Si of the mode i (i = a, b, c) is defined by |)3|] 

Si = -Trj(pjlnpi), (5.1) 

where pi is the reduced density operator of the i mode. In particular, the reduced density 
operator p a of the pump mode is 

p a = Tr 6c (p), (5.2) 

where p is the density operator of the whole system, and Tr 6c denotes tracing over the signal 
and idler variables. For the coherent state \y; k, L), we find 

Pa = jlidniK L)f\y\ 2n \n)(n\. (5.3) 

If S x and S y are the marginal entropies of two subsystems which compose the whole 
system with the total entropy S, the Araki-Lieb theorem can be expressed as 

\S x -S v \<S<S x + Sy. (5.4) 

In our case S — 0, and hence the marginal entropies of the pump and the signal-idler 
subsystems are equal: 
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S a = S bc . (5.5) 

Also, it follows from the Araki-Lieb theorem that 

\Sb - S c \ < S bc = S a < S b + S c . (5.6) 

A quantitative measure of the entanglement between two subsystems is the index of corre- 
lation 1513]: 



lory — S x -\- S y S xy . (5-7) 

In our case the index of correlation I a -bc between the pump and the signal-idler subsystems 
is equal to twice the marginal entropy of the pump mode: l a -bc = 2S a . 

Instead of evaluating entropies it is convenient to evaluate the purity parameters of the 
modes: 

Sr = 1-TW?). (5.8) 

It can be shown that the purity parameter Sf UT represents a lower bound for the correspond- 
ing entropy Si, i.e., Sf m < Si. Moreover, in the present case the purity parameters satisfy 
relations valid for the entropy, such as the Araki-Lieb theorem. In particular, this gives 

£P ur = o < S? ur = S^ r . (5.9) 

For the coherent state \y; k, L), we obtain 

sr = i - t-Jtxh EK(mfl!/l 4 " (5.io) 

I -VI • /V n=0 



This quantity is presented in Fig. [I] versus \z\ = l/\y\ for k = 1 and various values of L. 
As | z | increases, S% UT rapidly grows, reaches a maximum, and then decreases. The greater 
is the initial number L of the pump photons, the higher is the maximum of S% UT , i.e., the 
stronger is the entanglement. In the limit \z\ <C 1, the purity parameter increases as 

£P ur ^AkL\z\ 2 + 0{\z\' 1 ), (5.11) 

while in the limit \z\ ^> 1, it decreases as 

9 T 

We see that for \z\ ^> 1 the entanglement between the pump and the signal-idler subsystems 
is very weak. This result is expected due to the fact that the coherent states are obtained by 
the repeated application of the annihilation operator a with no creation operator a) being 
applied. Therefore, as \z\ becomes large, the pump mode tends toward the vacuum and is 
then more pure and more uncorrelated to the signal and idler modes. 

We also consider photon statistics of the pump mode a for the whole system in the 
coherent state \y; k, L). The photon-number distribution is given by 
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P a {n) = (n\p a \n) = ^ N 2 ^. (5.13) 

This distribution is plotted in Fig. ^ versus n for k — 1, L = 10, and various values of 
|z| = l/|y|. We see that P a (n) is very sensitive to the value of \z\. We can further investigate 
photon statistics of the pump mode by evaluating the expectation value and the variance of 
the number operator N a = a) a. It is easy to find the following expressions: 

<*•> = T A^< <» 4 > 

,, a a,j. r 2 d 2 A r dA ( r dA\ 2 

« AJV -) > = ^ + S*-Udrj ' (5 ' 15) 
where r = \y\ 2 and we have defined 

A(r; k, L) = rW(r; k, L) = Q^±_S$(_L ; 1 — L — 2fc;r). (5.16) 

We define 

n(r; k ,L) s ± in L; I — L — 2 k ; r) ] = L + ^ iTtl -"/-^r < 5 ' 17 ' 
Then the expectation value and the variance of iV a can be written as 

(4) = m, (5.i8) 

((AiV a ) 2 ) = [r 2 + (2A; + L)r]fi - r 2 fi 2 - Lr. (5.19) 

Using the integrals of motion, we can also infer the information about the signal-idler system: 

(K ) = k + L- (N a ), ((AK ) 2 ) = ((AN a ) 2 ). (5.20) 

Numerical results are presented in Fig. [| where (N a ) and ((AN a ) 2 ) are plotted versus 
\z\ = 1/ 'y/r for k = 1 and L = 10. We also consider some approximate results. In the limit 
\z\ ^ 1, we obtain 

(N a ) ^ L -2kL\z\ 2 + 0(\z\ 4 ), (5.21) 
((AN a ) 2 )^2kL\z\ 2 + 0(\z\ i ), (5.22) 

i.e., photon statistics of the pump mode is close to that of the Fock state \L). However, as 
\z\ increases, photon statistics rapidly changes. In the limit \z\ ^> 1, we obtain 

(N a ) « <(AiV a ) 2 > « L ^ Z 2 [ 2 _ 1 + C(kr 4 )- (5.23) 

Therefore, for \z\ ^> 1 photon statistics of the pump mode is close to be Poissonian. Because 
the purity parameter S% w in this limit is very small [see Eq. ( 5.12j) 1, i.e., the pump is almost 
disentangled from the signal-idler system, the state of the pump mode is close to a slightly 
excited pure state with Poissonian statistics. This result can be simply explained by the 
repeated application of the pump-mode annihilation operator with no creation operator 
being applied. For \z\ ^> 1 this leads to the strong depletion of the pump mode. 
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VI. LONG-TIME EVOLUTION AND FUNDAMENTAL LIMIT ON THE 

ENERGY TRANSFER 



We know that the coherent states \z;k,L) with \z\ <C 1 are produced from the initial 
state \L)\k, 0) during the first instants of the time evolution. However, this simple physical 
meaning is not valid in general for the coherent states with large values of \z\. Nevertheless, 
it should be interesting to study the connection between these states and the long-time 
evolution induced by the trilinear Hamiltonian. 

We consider the dynamics of the initial state \L)\k, 0). For general L the eigenvalue equa- 
tion ( |4.17| ) can be solved numerically by a diagonalization of the interaction Hamiltonian. 
For L = 1 and L = 2 this can be done analytically. By studying these simple situations, we 
can understand general features of the long-time dynamics. For L = 1 we obtain 

|V>t;M) =exv{-iH int t)\l)\k,0) = cos(Ait)|l)|fc,0) + sin(Ait)|0)|fc,l), (6.1) 



where k = i\K\, Ai = v2fc|/c|. We see that the evolution is periodical. Therefore, it is 
sufficient to consider only the time range < Ait < 2ir. The coherent state for L = 1 is 



, M> = |i)M) + V2MQ)IM) (62) 

/ l + 2k\z\ 2 



It is easy to check that this state can be written (up to an unimportant phase factor) in the 
form ( |6.1|) with real z = {2k)~ l l 2 tan(Ait). Thus the coherent state \z\ k, 1) exactly describes 
the dynamics induced by the trilinear Hamiltonian. The first instants of the time evolution 
(Ait <C 1) and the moments with Ait close to 7r and 2n are described by the coherent states 
with \z\ <C 1, while at the moments with Ait close to it/2 and 37r/2 the system evolves into 
the coherent state with \z\ 1. 

For L = 2 the situation is more complicated. In this case we obtain 

\^ t ;k,2) =exp(-zF int t)|2)|A;,0) = (4fc + l)- 1 { v /8fc(2fc + 1) sin 2 (A 2 t)|0)|/c, 2) 

+ y/2k(4k+ l)sin(2A 2 t) |1) | k, 1) + [1 + 4k cos 2 (A 2 t)] |2) | k, 0) }, (6.3) 



where A 2 = |\/8A; + 2\k\. The coherent state for L = 2 is given by 

2) = |2)|fc,0) + V4kz\l)\k, 1) + y/2k(2k + l)z 2 \0)\k, 2) 
^Jl + 4k\z\ 2 + 2k{2k+ l)|z| 4 

This state with \z\ 1 coincides with the state \i/)t',k,2) with A 2 t 1 (for z = |ac|£). 
However, for other moments of time this simple connection is not valid in general. For 
example, at the moment A 2 t = 7r/2 — e, £ <C 1, the state \ip t ; k, 2) of Eq. ( |6.3| ) has the form 



|V»t;A:,2> w (4A;+ l)- 1 [^8A;(2fc + 1)(1 -e 2 )|0)|fc,2) + ^/8fc(4fc + l)e|l)|fc, 1} 

+ (l + 4fe 2 )|2)|fc,0)l. (6.4) 
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On the other hand, the coherent state \z; k, 2) with \z\ ^> 1 can be written as 



\z; k, 2) « (1 - g 2 )\0)\k, 2) + V2g\l)\k, 1) + ^J(2k + l)/(2k)g 2 \2)\k, 0), (6.5) 



where we have assumed z to be real and positive, and g = l/(\ / 2k + lz). The state \z\ k, 2) 
of Eq. (|6.5| ) does not coincide in general with the state \ipt] k, 2) of Eq. ( |6.4| ). It is seen that 
these two states coincide only in the limit k ^> e~ 2 ^> 1 (with g = e). 

In order to understand why the system does not evolve in general into the coherent 
state \z;k, 2) with \z\ ^> 1, consider the result (|5.23 ) for the mean photon number (N a ). 



Indeed, for \z\ ^> 1, (N a ) 1S very small, i.e., almost all the photons are in the signal-idler 
system. This result contradicts to a fundamental limit on the energy transfer, that exists 
for the evolution of the initial state \L)\k,0). It was first shown by Drobny and Buzek 
| 56fl that such a fundamental limit on the energy transfer exists in the multiphoton down- 
conversion, when the signal mode b is initially in the vacuum state and the pump mode a is 
excited. A particular case is the two-photon down-conversion (i.e., the degenerate parametric 
amplification), described by the Hamiltonian ( [2.5|) with the two-photon SU(1,1) realization 
( [2.12| ). In this case only less than 3/4 of the pump photons can be transferred in average 
from the pump to the signal mode |56fl . 

A similar limit holds also for the evolution induced by the general Hamiltonian (|2.5|) for 



the initial state \L)\k, 0). For L — 1, we find 

(iV a )(t) = cos 2 (A 1 t). (6.6) 

In this case the single photon being initially in the pump mode can be completely transferred 
to the signal mode. Therefore, there is no prohibition on the production of the coherent 
states with arbitrarily large values of \z\. However, for L > 1 there does exist the limit on 
the energy transfer. For L = 2, we obtain 

7 8fcsin 2 (A 2 t) 8fcsin 4 (A 2 t) 

{Nam = 2 -^kTi (AkTTF- ( } 

The efficiency £ of the energy transfer is defined as the maximum mean number of photons 
which can be in principle transferred from the pump to the signal-idler system, divided by 
the initial number of the pump photons. For L = 2, the efficiency is 

In particular, in the degenerate parametric amplification we have k = 1/4 and £ = 3/4, 
in accordance with the results of Drobny and Buzek [Q. However, in the non-degenerate 
parametric amplification the efficiency £ rapidly approaches 1 as k increases. Therefore, the 
coherent states with large values of \z\ can be produced in the limit k 3> 1, as we have seen 
from the comparison of Eqs. ( |6.4| ) and ( |6.5| ). 

We also find an heuristic relation between the coherent states and the long-time evolution 
induced by the trilinear Hamiltonian. Consider the state \ipt] k, 2) at the moment A 2 t = vr/2: 

|2>|M> + ^/8fc(2fc+l)|0>|fc,2> 
\4>t;k,2) = 4kTl • (6 ' 9) 
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We compare this state with an even superposition of two coherent states, 



|2)|fc,0) + \ 2k(2k + l)z 2 \0)\k, 2) 
\z; k, 2) e = A(\z; k, 2) + | - z; k, 2)) = V , (6.10) 

y/l + 2k{2k + l)z A 

where A is a normalization factor and we assume z to be real. It is easy to see that the 
superposition state \z; k, 2) e coincides with the state \tp t ] k, 2) of Eq. (|6.9|) for z 2 = 2. In this 
special case the long-time dynamics can be described by using a discrete superposition of 
the coherent states. However, any quantum state obtained during the time evolution of the 
system can be expanded in the coherent-state basis by using the identity resolution (|4.1|). 



VII. CONCLUSIONS 



In this paper we have introduced the coherent states for quantum systems, whose dynam- 
ics is governed by a trilinear boson Hamiltonian. This dynamics is complicated because the 
operators which compose such a Hamiltonian do not form a finite-dimensional Lie algebra. 
Therefore, the standard formalism of Perelomov cannot be used in this case. Nevertheless, 
we have shown that in each invariant subspace there does exist an overcomplete set of the 
coherent states. During the first instants of the time evolution the system prepared in a 
special initial state will evolve into the coherent state. Using the coherent-state basis, we 
have constructed the resolution of the identity and the analytic representation in the com- 
plex plane. This analytic representation coincides with a double representation based on the 
Glauber coherent states of the pump and on the SU(1,1) Perelomov coherent states of the 
signal-idler system. We have studied the entanglement between the field modes and photon 
statistics of the coherent states. We have also considered connections between the coherent 
states and the long-time evolution induced by the trilinear boson Hamiltonian. 
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FIGURE CAPTIONS 



FIG. 1. The purity parameter SP ur of Eq. (fTIcD functions of \z\ = l/\y\ for the coherent 
states \y; k, L) with k = 1 and L = 1, 3, 10. The loss of the correlation for \z\ 3> 1 occurs because 
the pump mode tends in this limit toward the vacuum and is then more pure. 



FIG. 2. The photon-number distribution P a (n) of Eq. ( 5.13j ) for the coherent states \y; k, L) with 
k = 1, L = 10, and various values of \z\ = l/\y\. 



FIG. 3. The mean photon number (N a ) of Eq. ( |5.18 ) [curve (a)], and the variance ((AN a ) 2 } of 



Eq. (|5l9D [curve (b)] as functions of \z\ = l/\y\ for the coherent state \y;k,L) with k = 1 and 
L = 10. As \z\ increases, the pump mode is depleted and approaches the vacuum in the limit 
Izl — > oo. 
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